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ABSTRACT 

We study the physical and kinematic properties of the narrow line region (NLR) of the 
nearest obscured quasar MRK 477 (2^=0.037), using optical and near-infrared (NIR) 
spectroscopy. About 100 emission lines are identified in the optical+NIR spectrum (90 
in the optical), including several narrow optical Fe+ lines. To our knowledge, this is 
the first type 2 active galactic nucleus (AGN) with such a detection. The Fe+ lines 
can be explained as the natural emission from the NLR photoionized by the AGN. 
Goronal line emission can only be confirmed in the NIR spectrum. 

As in many other AGN, a significant correlation is found between the lines full 
width at half-maximum and the critical density log(ncrit)- We propose that it is caused 
by the outflow. This could be the case in other AGNs. 

The nuclear jet-induced ionized outflow has been kinematically isolated in many 
emission lines covering a broad range of ionization potentials and critical densities. It 
is concentrated within R ^fewxlOO pc from the central engine. The outflowing gas 
is denser (n >8 000 cm“^) than the ambient non-perturbed gas (n ^400-630 cm“^). 
This could be due to the compression effect of the jet induced shocks. Alternatively, 
we propose that the outflow has been triggered by the jet at R <220 pc (possibly at 
<30 pc) and we trace how the impact weakens as it propagates outwards following 
the radiation-pressure dominated density gradient. 

The different kinematic behaviour of [Fell] A1.644 /im suggests that its emission is 
enhanced by shocks induced by the nuclear outffow/jet and is preferentially emitted 
at a different, less reddened spatial location. 
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1 INTRODUCTION 

MRK 477 (SDSS J144038.1+533016, z=0.037 and luminos¬ 
ity distance Dl = 161 Mpe) is a type 2 luminous active galac¬ 
tic nucleus (AGN). Also k nown as I Zw 92, it was first 
identified bv IZwicl^ (|l966l ) as a compact galaxy. It is usu¬ 
ally referred to as the most l uminous Seyfert 2 in t he local 
Universe. As pointed out bv I Heckman et akl (|l997l ). it has 
the highest [O III]A5007 luminosity (L[oiii]=3.3x 10^^ erg 
s“^) of any of the 140 Seyfert nuclei (type 1 or 2) corn- 


emission lines - quasars: general - quasars: individual: 


piled by IWhittlel (|l992l ) and the fifth highest radio power 
(Fi. 4 GHz =2.0x10^° erg s“^ Hz“^ at 1.4 GHz). Indeed, its 
high L[oiii] places it in the regime of optically selected ob¬ 
scured qua sars (QS02), according to the selection criteria 
defined by IZakaniska et al.l (|200 L[oiii] >1.2x10^^ erg 
s~^. The [OIII] luminosity implies a bolom etric luminosity 
Lboi=S.6 xlO^^ erg s~^ (Stern & Laor |201^ : see also Lamas- 
tra e t al. l2009l ). which is in the quasar range (Shen et al. 
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The quasar host galaxy is interacting with an emission 
line companion, possibly a LINER, lo cated 50 arcsec (~36 
kpc) to the north (De Robertis et ah IiQSTIL Hubble Space 
Telescope (HST) images are shown in Fig.[JlJ. Both galaxies 
are connected by a faint bridge and prominent tidal tails 
are apparent. The ratio of stellar masses between the quas ar 
host and the companion is Ml/M2=1.6 (Koss et al. l2012l L 

The quasar host shows a compact blue central source, 
peculiar among type 2 AGNs, but no evidence for a broad 
line region (BLR), prompting some authors to suggest that 
the quasar lacks one. However, polarization observ ations 
have revealed a hidd en BLR (e.g. Tran et al. Il992l : Tran 
I 1995 I : Shu et al. l2007l L The nuclear continuum polarization 
is <1.5%, significantly lower than that of the broad lines, 
implying that the majority of the optical continuum cannot 
be due to scattered light from the hidden type 1 nucleus. 
The blue central source is instead associated with a compact 
dusty starburst (effective radius<0.2 kpc) which occurred 
~6 Myr ago in a brief period of time and whose light domi¬ 
nates fro m the ultraviolet to the near-i nfrare d (NIR; Heck¬ 
man et al. ll997l . Gonzalez Delgado et al. ll998l L This is there¬ 
fore a starburst-AGN hybrid system. MRK 477 has a star¬ 
forming rate SFR=24 M© yr“^ inferred from the infrared 
(IR) luminosity (assumin g tha t it is starburst-dominated) 
log ^^ = 11.14 (Krug et al. l2010l L This places it in the regime 

of luminous IR galaxies (LIRGs: 11^ log'll <12). The de¬ 
tection of the Wolf Rayet (WR ) blue bump around the He 
HA4686 line bv I Heckman et all (|l997l j led these authors to 
propose that MRK 477 is a luminous member of the class of 
WR galaxies. 

The QS02 host harbors a supermassive black hole with 
a mass in the range log(MBH)~7.18 (as derived from the 
Mbh versus cr* correlation) to 8.84 (as in ferred from the 
polarimetric H/3 line; Zhang, Bian & Huang [2008l j . 

Although radio quiet according to the L[oiii) versus 
Trad, radio loud versus radio ciuiet, classification criteria (ex¬ 
plained in Villar-Martm et al. l2014l L MRK 477 shows a clear 
excess of radio emission compared to that expected from the 
stellar contribution. The 8.4 GHz Very Large Array (VLA) 
radio continuum map (0.26 arcsec resolution) shows a triple 
radio source (~1.2 arcsec total extension) whose morphol¬ 
ogy correlates with that of the narrow line region (NLR). 
The size and overall northeast-southwest axis of the radio 
structures are shared by the [OHI] emission, which extends 
up to a similar distance to the north-east (Heckman et al. 

IShuder fc OsterbrockI (|l98lh reported for the first time 
that the emission lines have a large blue-ward asymmetry 
(i.e., excess) at low i ntens ity level relative to the continuum. 

I Villar Martin et m] (|2014l i proposed that the ionized nuclear 
outflow responsible for such an asymmetry has been trig¬ 
gered by the interaction between the radio source and the 
NLR. These authors proposed that negative feedback can be 
triggered by the radio structures in a significant fraction of 
radio quiet quasars. Thanks to its high intrinsic luminosity 

^ The HST images were retrieved from the HST Science Legacy 
Archive. They are part of the programmes 9379 with principal 
investigator PI H. Schmitt (ACS/HRC image), 8597 with PI M. 
Michael Regan (WFPC image) and 7330 with PI J. Mnlchaey 
(NICMOS image). 


and closeness, MRK 477 is an excellent test object to gain 
further insight into this mechanism. 

We present here a detailed study of the optical and 
near-infrared spectra of MRK 477. We explore a diversity 
of aspects that provide a more complete understanding of 
the nature of this object, example of a type 2 quasar in the 
nearby universe, as well as a starburst-AGN hybrid system. 
The paper is organized as follows. The spectra are described 
in Sect. 2, together with the spectral fitting procedure. We 
present results in Sect. 3 regarding line identification, cor¬ 
relations between the gas kinematics and a diversity of pa¬ 
rameters (ionization potential, critical density), the spatial 
extension of the ionized gas, the ionized outflow, the pres¬ 
ence of narrow Fe^ emission lines, the coronal line spectrum 
and the detection of WR features. The results are discussed 
in Sect. 4 and the conclusions are presented in Sect. 5. 

We assume Qa=0.7, Qm= 0.3, iLo=71 km s“^ Mpc~^ . 
At z =0.037, 1 arcsec corresponds 0.725 kpc fWright [20061 ) . 

2 DATA SET 

2.1 Optical spectrum 

The optical spectrum was obtained as par t of the Sloan Dig¬ 
ital Sky Survey (SDSS, York et al. I 2 OOOI L It spans the rest 
frame range ~3660-8880 A range (Figs. [2] and [3)). It cor¬ 
responds to an aperture defined by the 3 arcsec diameter 
SDSS hbre (~2.2 kpc at z =0.037) centred at the galaxy 
nu cleus. For comparison , the triple radio source identified 
bv iHeckman et J (|l997l i has a total extension of ~1.2 arc¬ 
sec and thus it is well within the fibre area. 

The spectral resolution is ~180±20 km s“^. All values 
of emission line full width at half-maximum (FWHM) have 
been corrected for instrumental broadening. 

The main emission lines have large equivalent widths 
and the stellar features are comparatively weak. Underlying 
stellar absorption of the Balmer lines is expected to be neg¬ 
ligible and therefore subtracting the stellar continuum is not 
necessary for our purposes. 

2.2 NIR spectrum: observations and data 
reduction 

NIR H+K long-slit spectra were obtained with the NIR cam¬ 
era/spectrometer LIRIS (Long-slit Interm ediate Resolution 
Infrar ed Sp ectrograph; Acosta-Pulido et al. l2003l : Manchado 
et al. l2004| j . attached to the Gassegrain focus of the 4.2 
m William Herschel Telescope (WHT). LIRIS is equipped 
with a Rockwell Hawaii 1024x1024 HgGdTe array detector, 
whose spatial scale is 0.25 arcsec pixeU^. 

The NIR observations were performed as part of 
the Isaac Newton Group service programme (programme 
SW2014b22) in two different nights, 28 March 2015 and 02 
May 2015. The H+K grism was used. It covers the 1.388- 
2.419 jJLVsi spectral range and provides a dispersion of 9.7A 
pix“^. The seeing size during the March observations was 
FWHM~2"(as measured from the standard star). A 1" wide 
slit was used oriented at position angle PA=43° East of 
North. This PA is aligned with the radio ax is and the NLR 
axis a s seeing in the HST image presented bv iHeckman et all 
(Il997h . The resulting spectral resolution was 35.7+2.2 A (or 
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Figure 1. Hubble Space Telescope (HST) images of MRK 477. They were obtained with Wide-Field Planetary Camera 2 (WFPC2, 
F606W filter, left), the Advanced Camera for Surveys (ACS/HRC, F330W filter, middle) and the near-infrared Camera and Multi-Object 
Spectrometer (NICMOS, F160W filter, right) respectively. The flux intensity maps are represented in logarithmic scale such that F [erg 
s“^ cm“^ arcsec”^] = where A is a constant factor (i.e., -18.3, -14.5 and -16.4, respectively, for WFPC2, ACS and NICMOS 

images) and b is indicated in the coloured bar. The white box in the WFPC2 image identifies the FoV covered by the ACS and NICMOS 
images. The horizontal line at the bottom left of these images corresponds to a scale of 10 kpc for the WFPC2 image and to a scale of 2 
kpc for the ACS and NICMOS images. The spatial scale for this galaxy is 0.725 kpc/ arcsec . North is at the top and East to the left in 
all the panels. 


550zb34 km at the observed wavelength of Paa). In Fig.[4| 
we show the NIR spectrum extracted from a 3.5" aperture, 
centred at the position of the nucleus. The total exposure 
time on source was 6400 s (16x400 s). 

The May data were obtained with seeing 
FWHM=0.65zb0.05"(measured from the standard star). 
A 0.75" slit was used, also with PA=43°. The resulting 
spectral resolution was 23.2zbl.2 A (or 318zbl8 km s“^). 
The observations were performed through clouds and the 
signal to noise ratio of the final spectrum was low. In 
spite of this, the good seeing conditions allowed to obtain 
more accurate information on the spatial extension of the 
strongest emission line Paa. The total exposure time on 
source was 4000 s (10x400 s). All FWHM values of the NIR 
spectra have been corrected for instrumental broadening. 

The data were reduced following standard procedures 
for near-IR spectroscope using the LIRISDR dedicated soft¬ 
ware within the IRARj enviromen t. For a detailed de¬ 
script ion of the reduction process see iRamos Almeida et al.l 
(|2009l ). Consecutive pairs of AB two-dimensional spectra 
were subtracted to remove the sky background. The result¬ 
ing frames were then wavelength calibrated and flat-fielded 
before registering and co-adding all frames to provide the 
final spectra. 

The absolute flux calibration of the March spectrum is 
intended to be an approximation since the spectra of the 
comparison star are likely subjected to slit losses due to 
centring and tracking errors. We note that in this work we 
refer to flux ratios only, and therefore the uncertainty in the 


^ IRAF is distributed by the National Optical Astronomy Obser¬ 
vatory, which is operated by the Association of Universities for 
the Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation (http://iraf.noao.edu/). 


absolute flux calibration of the spectra does not affect our 
results. 


2.3 Spectral fitting 

In order to study the kinematic and physical properties of 
the nuclear ionized outflow we have fitted the spectral pro¬ 
files of numerous emission lines (Sect. 3). For this we used 
the STARLINK package DIPSO. This software is based on 
the optimization of fit coefficients, in the sense of minimizing 
the sum of the squares of the deviations of the fit from the 
spectrum data. The output from a completed fit consists of 
the optimized parameters (Gaussian central A, FWHM, in¬ 
tensity peak, flux) and their errors (calculated in the linear 
approximation, from the error matrix). 

Two methods were attempted. In method I, we assumed 
that all emission lines have the same kinematic substructure 
as [OIII]AA4959,5007 (same number of kinematic compo¬ 
nents with identical FWHM in km s“^ a nd rel ative veloc¬ 
ities, as explained in Villar Martin et al. l2014l ). In physi¬ 
cal terms, this method corresponds to a situation where all 
gaseous regions emit all lines, although with different rela¬ 
tive fluxes. 

On the other hand, this is not necessarily the case and 
some lines might come from different regions, possibly re¬ 
sulting in variations of the kinematic substructure from line 
to line. As an example, low critical density lines (~fewxl0^ 
cm“^) are quenched in high density gas with n >10® cm“^. 
To account for this possibility, whenever possible we have 
also attempted to fit the lines without applying prior kine¬ 
matic restrictions (method H). 

Also, all fits that produced unphysical results were re¬ 
jected even if mathematically valid (e.g. a fit producing [NH] 
very different from the theoretical value 3.0). 

As we will see the lines are complex with multiple kine- 
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Figure 2. SDSS spectrum of MRK 477. Flux is given in units of xlO erg s ^ cm ^ A and the rest-frame wavelength is in A. 
Some of the strongest emission lines are indicated. 






Figure 3. Four zoomed spectral windows covering the full spectral range are shown to highlight some of the weakest lines. Unless 
otherwise specified, fiuxes in these and other figures are given in units of erg s“^ cm“^ A“^. Some relevant lines are indicated. 
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matic components and they are often severely blended with 
neighbour lines so that the hts are sometimes complicated 
and not unambiguous. Applying methods I and II whenever 
possible we account for the uncertainties involved in a more 
realistic way than using a single method. 

Both methods do not always produce accept¬ 
able results. For instance, the [SII]AA4068,4076+H(5 and 
[SII]AA6716,6731 blends could not be successfully htted ap¬ 
plying full constraints from the [OIII] lines. The opposite 
occurs with the Ha+[NII] blend. 

Methods I and II could be successfully applied to Hy, 
[OIII]A4363, H/3 and [OI]A6300. It is found that, [OIII]A4363 
is the most uncertain. While the narrow and intermediate 
components isolated in the hts (see below) are consistent 
within 30 and 10 per cent, respectively, with the two meth¬ 
ods, a discrepancy of a factor of ~60 per cent is found for 
the broadest component. The results for all three kinematic 
components isolated in the hts of the other three lines dif¬ 
fer by <20 per cent. These uncertainties will be taken into 
account when relevant. 


3 RESULTS 

3.1 Line identification. 

We identify ~90 emission lines in the SDSS optical spec¬ 
trum of MRK 477 (Table many of them detected for 
the hrst time in this object. For some features the iden- 
tihcation is ambiguous and/or it could be the contribu¬ 
tion of several emission lines. In such cases, all the pos¬ 
sible identihcations or contributors are quoted separated 
with As an example, a broad feature with an asym¬ 
metric prohle is detected at A ~5272 A. This is identihed as 
“[Fe III];[Fell];[Fell];[Fe VII]” because it might be the blend 
of [FeIII]A5270.4, [FeII]AA5273.4,5276.0 and [FeVII]A5276.4. 
When two or more lines are known to contribute to a given 
feature, all are quoted and separated with • For instance, 
the [Oil] at A ~3727 and [NI] at A ~5200 doublets are shown 
as “3727.0+3728.8” and “5197.9+5200.4”. 

The works by IVeron-Cettv. Jolv fc VeronI (|2004l j and 
IVeron-Cettv et all (|2013l l have been used for iron (Fe) line 
identihcations. The critical densities Ucrit of forbidden tran¬ 
sitions are quoted in column 5 of Table [T] when avail¬ 
able. Most have been computed with the PYNEB pack¬ 
age optimized for the analy sis of emission lines (Luridi- 
ana, Morisset V Shaw l2014l b These have been estimated 
for Te ^15,000 K. The rest h ave been retrieved from 
IPe Robertis fc Osterbroclj (|l984l b 

The lines identihed in the NIR spectrum are shown in 
Table O 


3.2 Relationships between the kinematics and 

T^crit ? IPlow and IPhigh* 

We plot in Fig. [5] the FWHM (top panels) and the velocity 
shift Vs (bottom) versus the critical density {ricrit, Table [J), 
the lower (IPiow) and upper (IPhigh) ionization potentials for 
a subsample of ~28 optical and 2 NIR forbidden emission 
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MRK477 NIR spectrum 



Figure 4. NIR LIRIS-WHT spectrum of MRK 477 extracted 
from a ~l"x3.5'' aperture centered on the nucleus. The detected 
lines are indicated. Flux is in units of erg s“^ cm“^ 

and A in /rm. 

lines 0 Vs corresponds to the difference between Aair and Aobs 
in Table [H where Aair is the vacuum wavelength and Aobs 
is that measured from the spectrum corrected for z. Only 
lines for which both parameters could be measured with rea¬ 
sonable accuracy are plotted. As an example, this was not 
the case for [SiX]A1.430 and [Fell] A 1.257 due to the weak¬ 
ness of the lines and/or the distortion of the line profiles by 
sky residuals. The optical and NIR lines span a wide range 
in critical densities (~1.3x 10^-6.3x 10^ cm“^) and ioniza¬ 
tion potentials (IPhigh ~13-205 eV). A significant correla¬ 
tion is found between the FWHM and the critical density 
(Spearman correlation coefficients rs=0.79 and p=0.000001, 
excluding [Fell]Al.644), as already found by De Robertis V 
Osterbrock (1986) whose data are shown as open circles. No 
signihcant correlation is found with the ionization potentials 
(rs=0.43 and p=0.02 for FWHM versus IPiow and rs=0.42 
and p=0.03 for FWHM versus IPhigh)- No trends are found 
either for Vs except maybe an apparent preference for red- 
shifts (Vs >0) for lines with the highest critical densities. 

The NIR [Fell]Al.644 line is a clear outlier in Fig. [5l 
being too broad compared to what is expected for its critical 
density (5.6x10^ cm“^) and ionization potentials (7.9 and 

16.2 eV) according to the general trends dehned by the other 
lines. We will discuss this in more detail below. 

The correlation between FWHM and Ucrit has been 
found in many type 1 and type 2 AGNs and it holds for 
AGN l umino sities that differ by a factor of up to ~5000 
(Espev H^ . Es, on the other hand, shows a different be- 
havio r from object to object (e.g. Appenzeller V Ostreicher 
Il988h . 

The FWHM versus Ucrit correlation implies that the 
line emission from the NLR originates in different subregions 
under different physical conditions and kinematic proper¬ 
ties with a broad density range. Higher critical density lines 
(many of which have also very high ionization potentials) are 
broader because they are emitted predominantly by higher 

^ The FWHM and Vs values have been derived from single Gaus¬ 
sian fits to the emission lines. 
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(1) 

Species 

(2) 

^air 

(A) 

(3) 

^obs 

(A) 

(4) 

Flux 

(5) 

(cm-3) 

Flux(H/3) 

[Oil] 

3727.0 + 3728.8 

3726.9 

2.26 + 0.18 

(1.3/4.5) X 10^ 

H12 

3750.2 

3750.1 

0.014 + 0.003 


[FeVII] 

3758.9 

3759.3 

0.046 + 0.005 

3.7X lo"^® 

Hll 

3770.6 

3770.4 

0.023 + 0.002 


HIO 

3797.9 

3797.5 

0.026 + 0.003 


H9 

3835.4 

3835.5 

0.045 + 0.003 


[Nelll] 

3868.8 

3868.6 

1.04 + 0.03 

1.3X lo"^ 

Hel + H8 

3888.7+3889.05 

3888.8 

0.16 + 0.02 


[NeIII]+H7 

3967.5 + 3970.1 

3968.2 

0.43 + 0.04 

1.3X lo"^ 

Hel 

4026.2 

4026.1 

0.013 + 0.002 


[SII] 

4068.6 

4068.6 

0.19 + 0.01 

3.Ox 10® 

[SII] 

4076.4 

4076.1 

0.065 + 0.03 

1.5X 10® 

H5 

4101.7 

4101.5 

0.24 + 0.02 


[FeV] 

4180.9 

4180.7 

0.008 + 0.002 


FeII];[FeV] 

4227.2;4228.0 

4228.2 + 0.1 

0.017 + 0.002 


[Fell] 

4244.0 

4244.5 

0.019 + 0.002 


[Fell] 

4276.8 

4278.3 + 0.1 

0.013 + 0.002 


[Fell] 

4287.4 

4289.1 + 0.4 

0.014 + 0.002 


H7 

4340.5 

4340.31 

0.45 + 0.05 


[OIII] 

4363.2 

4363.2 

0.23 + 0.02 

3.4X lo"^ 

[Fell] 

4416.3 

4416.4 

0.021 + 0.002 


Hel 

4471.4 

4472.2 

0.035 + 0.002 


[Fell] 

4570.0 

4570.9 

0.018 + 0.005 


[FeIII];CIV 

4658.1 

4658.7 

0.034 + 0.001 


Hell 

4685.7 

4685.7 

0.138 + 0.008 


[Felll] 

4701.5 

4700.2 

0.010 + 0.008 


[ArlV] 

4711.4 

4711.5 

0.011 + 0.002 

1.7X 10^ 

[NelV] 

4725.5;4726.8 

4725.1 

0.008 + 0.003 

(6.2/2.5) X lo"^ 

[ArlV] 

4740.2 

4740.2 

0.016 + 0.003 

1.6X 10® 

[Felll] 

4754.7 

4754.0 + 0.5 

0.007 + 0.002 


[Felll] 

4769.5 

4769.6 

0.006 + 0.002 


[Felll] 

4777.7 

4777.5 

0.007 + 0.002 


[Fe II] 

4814.5 

4816.3 + 

0.016 + 0.002 


H/3 

4861.3 

4861.3 

1.00 


[FeVII] 

4893.4 

4893.9 

0.008 + 0.002 


[OIII] 

4958.9 

4958.9 

3.7 + 0.3 

8.6X 10® 

[OIII] 

5006.8 

5006.8 

10.5 + 0.3 

8.6X 10® 

[FeVI] 

5145.8 

5145.9 

0.007 + 0.002 

2.7X lo"^® 

[FeVII] 

5158.9 

5159.4 

0.042 + 0.005 

3.3X 10®“ 

[FeVI] 

5176.4 

5177.0 

0.010 + 0.001 

3.3X 10"^“ 

[NI] 

5197.9 + 5200.4 

5199.2 

0.05 + 0.01 

(1.6/0.51) X 10® 

[Fe III]; [Fell]; [Fell]; [Fe VII] 

5270.4; 5273.4; 5276.0; 5276.4 

5272.1 

0.058 + 0.004 


[CaV] 

5309.1 

5309.7 

0.032 + 0.002 

5.1X10’^ 

[FeVI] 

5335.2 

5335.5 

0.017 + 0.002 


[Fell] 

5379.0 

5379.7 

0.011 + 0.002 


Hell; [Fell] 

5411.5;5412.7 

5412.0 

0.015 + 0.005 


Fell]; [FeVI] 

5425.3;5425.7 

5424.1 

0.018 + 0.006 


[FeVI] 

5485.0 

5484.5 

0.005 + 0.001 


[Cl III] 

5517.7 

5517.5 

0.006 + 0.002 

8.5X 10® 

[ArX];[FeII];FeII] 

5533.0;5527.6;5534.8 

5531.39 

0.013 + 0.003 


[Cl III] 

5537.9 

5537.9 

0.009 + 0.002 

2.9X 104 

[OI] 

5577.34 

5577.9 

0.014 + 0.001 

9.0X lo"^ 

[FeVI] 

5631.1 

5630.9 

0.008 + 0.002 


[FeVI] 

5677.0 

5677.1 

0.005 + 0.001 


[FeVII] 

5720.9 

5721.8 + 0.2 

0.047 + 0.002 

3.6X lo"^® 

[Nil] 

5754.6 

5754.7 

0.040 + 0.002 

1.8X lo"^ 

Hel 

5875.6 

5875.7 

0.131 + 0.007 


[FeVII] 

6086.9 

6087.7 

0.12 + 0.01 

3.6X 10"^“ 

[OI] 

6300.3 

6300.5 

0.81 + 0.02 

1.1 X 10® 

[Sill] 

6312.1 

6311.8 

0.016 + 0.003 

1.6X lo"^ 

[OI] 

6363.8 

6365.0 

0.26 + 0.02 

1.1 X 10® 

[ArV] 

6435.1 

6435.7 

0.010 + 0.002 

1.4X lo"^ 

[Nil] 

6548.1 

6548.1 

0.53 + 0.05 

1.2X 10® 

Ha 

6562.8 

6562.8 

3.9 + 0.2 


[Nil] 

6583.5 

6583.4 

1.5 + 0.1 

1.2X 10® 

Hel 

6678.2 

6678.7 

0.035 + 0.005 


[SII] 

6716.4 

6716.5 

0.60 + 0.03 

1.5X 10® 

[SII] 

6730.8 

6730.9 

0.75 + 0.03 

4.Ox 10® 

[ArV] 

7005.8 

7005.7 

0.014 + 0.003 

1.3X lo"^ 

Hel 

7065.7 

7065.3 

0.07 + 0.005 


[Arlll] 

7135.8 

7135.8 

0.20 + 0.01 

6.4X 10® 

[Fell] 

7155.2 

7156.8 

0.06 + 0.01 


[Fell] 

7172 .0 

7173.8 

0.014 + 0.004 


[Call] 

7291.5 

7292.8 

0.02 + 0.01 


[Oil] 

7318.9 

7321.0 

0.23 + 0.04 

6.5X 10® 

[Oil] 

7339.7 

7329.2 

0.32 + 0.04 

6.5X 10® 

[Nill] 

7377.8 

7379.1 

0.040 + 0.008 


[Fell] 

7452.6 

7453.9 + 0.5 

0.016 + 0.002 


FeI;[FeII] 

7634.0;7637.5 

7639 + 2 

0.007 + 0.002 


[Arlll] 

7751.1 

7750.9 

0.062 + 0.003 

6.4X 10® 

Mgll; [Nilll]; [FeXI]; Mgll 

7877.0; 7889.9; 7891.8;7896.0 

7888.8 

0.023 + 0.003 


Mgll 

8232.0 

8232.7 + 0.5 

0.020 + 0.005 


OI 

8446.0 

8446.6 

0.029 + 0.002 


HI Pa 14 

8598.5 

8595.4 

0.010 + 0.003 


[Fell] 

8605.5 

8604.5 

0.033 + 0.004 


[Fell] 

8617.0 

8617.8 

0.05 + 0.02 


HI Pa 12 

8750.6 

8750.2 

0.013 + 0.002 



Table 1. Emission lines identified in the SDSS spectrum of MRK 477. The observed line wavelengths (column 3) have been 

determined by fitting a single Gaussian to the emission lines. Errors are quoted when they provide relevant information regarding 
possible multiple identifications or there is a significant shift between the observed and the air wavelengths (column 2). The line ratios 
(column 4) are given relative to H/3, with F(H/3)=(1.01±0.04)x 10“^^ erg s“^ cm“^. They are not corrected for reddening. The critical 
densities are quoted when available. Superscript “ in column (5) indicates that ricrit has been taken from De Robertis & Osterbrock 
(1984). The other ricrit values have been computed with the PYNEB software (see the text). 
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(1) 

Species 

(2) 

^air 

(A) 

(3) 

-^obs 

(A) 

(4) 

Flux 

(5) 

Tl'crit 

(cm-3) 

Flux(Br'y) 

[SiX] 

1.4305 

1.4308 

0.61T0.10 

6.3xl0®(“) 

[Fe II] 

1.5339 

1.5337 

0.79T0.10 

4.6x104('’) 

[Fe II] 

1.6440 

1.6435 

4.5T0.5 

5.6x104('’) 

PaQ 

1.8756 

1.8746 

12.6T1.4 


Ha 1-0S(3) 

1.9576 

1.9564 

0.84T0.10 


[Si VI] 

1.9629 

1.9626 

0.97T0.12 

6.3xl0®(“) 

Ha 1-0S(2) 

2.0338 

2.0329 

0.20T0.02 


Hel 

2.0587 

2.0575 

0.29T0.03 


Ha l-OS(l) 

2.1218 

2.1206 

0.63T0.17 


Br7 

2.1661 

2.1650 

1.00T0.15 



Table 2. Emission lines identified in the NIR spectrum of of MRK 477 extracted from a arcsec X 3.5 arcsec aperture centered on 
the nucleus, ^^^ricrit from Rodriguez Ardila et al. (2011). ^^^ricrit from Drougados et al. (2010). Fluxes are quoted relative to Br 7 , for 
which F(Br7)=(3.0±0.3)X10“^^ erg s“^ cm“^, although the absolute flux calibration is uncertain (see the text). 
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Figure 5. The FWHM (top panels) and the velocity shift Vs (bottom panels) of a subsample of ~30 forbidden lines are plotted against 
the critical de nsity ( log(nr,rU,)). the lower and higher ionization potentials of the species involved. The measurements from De Robertis 
& Osterbrock (|l986l . dR & Ost 86) are also shown in the top-left panel as open circles. The outlier at log(ncrit) ~3.1 and FWHM~500 
km s“^ corresponds to [OIIJA3727, which appears artificially broader due to blended doublet components. The NIR lines are marked 
with red triangles. [Fell] Al.644 is a clear outlier in the top panels. It is too broad compared to what is expected from the general trend 
defined by the rest of the lines. 
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8 Villar Martin et al. 


Comp. 

FWHM 

(km s“^) 

(km s“^) 

Fi 

F(H/3) 

[SIIJA4068 

Narrow 

Interm. 

Broad 

94±37 

505±35 

1534±97 

-2±8 

17±10 

-376±85 

0.07±0.01 

0.26±0.02 

0.25±0.05 

[SIIJA4076 

Narrow 

Interm. 

Broad 

94±37 

505±35 

1534±97 

-2±8 

17±10 

-376±85 

0.010±0.006 

0 .12±0.02 

0.08±0.07 

H<5 

Narrow 

Interm. 

Broad 

94±37 

505±35 

1534±97 

-2±8 

17±10 

-376±85 

0.23±0.04 

0.23±0.06 

0.26±0.03 

H7 

Narrow 

Interm. 

Broad 

157±28 

529±40 

1830±133 

-7±7 

36±10 

-224±51 

0.47±0.03 

0.38±0.04 

0.57±0.06 

[OIII] A4363 

Narrow 

Interm. 

Broad 

157 ±28 
529±40 
1830±133 

-7±7 

36±10 

-224±51 

0 .11±0.02 

0.43±0.03 

0.27±0.08 

H/3 

Narrow 

Interm. 

Broad 

137±22 

471±22 

1527±75 

-2±6 

7±7 

-188±24 

(4.1±0.4) X 
(3.7±0.1) X 10“^^ 
(2.3X0.1)X10“1^ 

[OIII]A5007 

Narrow 

Interm. 

Broad 

177±15 

545±15 

1839±53 

5±9 

22±11 

-227±23 

8.5±0.9 

13.1±0.5 

9.6±0.7 

[OI]A6300 

Narrow 

Interm. 

Broad 

211±11 

563±29 

1745±40 

-4±5 

17±6 

-257±24 

0.60±0.07 

0.91±0.08 

1.08±0.07 

Hcc 

Narrow 

Interm. 

Broad 

As [OIII] 

As [OIII] 

As [OIII] 

As [OIII] 

As [OIII] 

As [OIII] 

3.1±0.3 

4.9±0.2 

4.2±0.3 

[NII]A6583 

Narrow 

Interm. 

Broad 

As [OIII] 

As [Olllj 

As [OIII] 

As [OIII] 

As [Olllj 

As [Olllj 

1 .0±0.1 

1.70±0.07 

0.9±0.2 

[SII]A6716 

Broad 

Interm. 

Broad 

186±9 

501±17 

1397±76 

-3±5 

-4±6 

-493±89 

0.67±0.07 

0.68±0.04 

0.30±0.08 

[SII]A6731 

Narrow 

Interm. 

Broad 

186±9 

501±17 

1397±76 

-3±5 

-4±6 

-493±89 

0.66±0.07 

0.80±0.03 

0.77±0.08 


Table 3. Results of the spectral decomposition of the main lines 
in the spectrum of MRK 477. Three kinematic components are 
isolated in all cases, including a broad blueshifted component 
emitted by the ionized outflow. No values are quoted for Ha and 
[Nil] because successful fits could only be obtained with full con¬ 
straints from the [OIII] lines (see Sect. 2.1). The flux of each com¬ 
ponent relative to H/3 is given in the last column. The exception 
is H/3, for which the actual fluxes are quoted. 


velocity, high density and high ionization regions closer to 
the central engine. 

We have done n estimations using the density- 
diagnostic line ratios [SII]A6716/A6731, [Cl III]A5517/A5537 
and [ArIV]A4711/A4740, which are not affected by ex¬ 
tinction or Te variations. They are sensitive to different 
n regimes. The [SII] ratio (ncrit=l-5x 10^ and 4.0x10^ 
cm“^, respectively) is sensitive to n variations ~severalx 10^ 
cm“^. These lines are suppressed in high density re- 



Shift (arcsec) 

Figure 6. Spatial profiles (normalized fluxes) of the NIR contin¬ 
uum and Paa along PA 43° N to E, compared with the seeing 
disk (star). Both the continuum and Paa are spatially extended. 


gions by collisional de-excitation. Other line pairs that 
can be overcome this limitation are [Cl III] A5517/A5537 
(ricrit—S.bxlO^ and 2.9x10^ cm“^), sensitive at n 
cm“^ and [ArIV]A4711/A4740 (ricrit— 1.7x10^ and 1.6x10^ 
cm “^) at n ~10^-10^ cm“^. 

For MRK 477, [SII]A6716/A6731=0.70±0.06, 

[C1III]A5517/A5537=0.74±0.22 and [ArIV]A4711/A4740 
=0.67zb0.09, impling densities n=1945l4gQ cm“^ (from 
[SII]), 725015100°° (from [Cl III]) and n=14,700l375o (from 
[ArlV]). So, the [SII] doublet implies the existence of gas 
with n ~2000 cm“°, the [Cl III] doublet, in spite of the 
larger uncertainties suggests higher densities and the [ArlV] 
doublet confirms the existence of gas with n >10^ cm“^ in 
the NLR of MRK 477. 


3.3 Spatial extension of the ionized gas 

The only available information about the spatial distribu¬ 
tion of the ionized gas in MRK 477 is prov ided by the [OIII] 
image presented bv lHeckman et aP (|l997l ) . It shows a bright 
knot or ridge of emission at ~0.4 arcsec to the northeast of 
the central source and aligned with the radio axis. This knot 
is not resolved from the nucleus in our NIR spectrum, for 
which the slit was roughly aligned along the same axis. We 
find, on the other hand, that Paa is extended at both sides 
of the central source. We show in Fig. [6] the spatial pro¬ 
files of (1) a star observed during the May run (seeing size 
FWHM=0.65±0.05), (2) the MRK 477 continuum and (3) 
the Pa flux, with the underlying continuum subtracted. Al¬ 
though more compact than the continuum, Paa is extended, 
as demonstrated by the clear excess above the seeing wings, 
up to a maximum radial distance of ~2 arcsec or ~1.5 kpc 
to the NE. 


3.4 The ionized outflow 

IVillar Martin et"al] (|2014l ) performed a kinematic and ion¬ 
ization characterization of the nuclear ionized outflow 
in MRK 477. Its signature is a broad blueshifted kine¬ 
matic component in the [OIII]AA4959,5007 lines, with 


© 2002 RAS, MNRAS 000 , [I]-?? 






































The NLR of the nearest obscured quasar 9 




Figure 7. Left: data (black), fits (cyan) and residuals (dashed-red) for a diversity of emission lines. Right: individual kinematic compo¬ 
nents isolated in each line. Different line styles are used for different kinematic components. The same red colour is used for components 
with the same redshift: long-dashed red (narrow component), dotted-red (intermediate component), dot-dashed blue (broad component; 
i.e. the outflow emission). From top to bottom: [SII] AA4068,4076 and H(5; Hy and [OIII] A4363; [Nil] AA6548,6583 and Ha; [SII] AA6716,6731. 
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10 Villar Martin et al. 



Figure 8 . e = parametrizes the relative contribution between the outflow flux emission and the rest of the line flux. The 

correlation with the critical density (panel A) implies that the outflow emission is relatively stronger in higher density gas. [OI] ‘a’ and 
‘b’ correspond to the A6300 and A6312 lines respectively; [NeIII]‘a’ and ‘b’ to A3869 and A3967; [SII] ‘b’ and ‘r’ refer to the blue and red 
[SII] doublets at AA6048,6076 and AA6716,6731 respectively. [Fell]A1.644p,m (red triangle) is again an outlier. 


FWHM~1850 km s~^ and Vs ~-225 km s~^ . Based on the 
high turbulence of this gaseous component, the line ratios 
consistent with AGN (rather than stellar) excitation pro¬ 
cesses and the relatively high contribution of the outflowing 
gas to the total line fluxes, the authors proposed that it has 
been triggered by the 1.2 arcsec scale radio source. This is 
supported by the correlation between the radio and [OIII] 
morphologies (Heckman et al. 1997), which demonstrates 
that the radio source is interacting with the NLR. 


5.^.7 Isolating the outflow emission in a diversity of 
emission lines 

By means of the kinematic decomposition of the spectral 
profiles, we have isolated the emission from the outflowing 
gas in numerous emission lines (Sect. 2.1). Our goal is to 
characterize its physical properties and to constrain more 
accurately its spatial location. 

We show in Table [3] the results of the fits for several 
optical emission lines (see also Fig. [7]). FWHM and Vs cor¬ 
respond to those fits where no prior constraints from [OIII] 
were applied (method II, see Sect 2.3). When both methods 
could be applied, the flux values and errors account for the 
dispersion allowed by them. 

Generally, we find rather coherent results for all strong 
lines for which multiple component fitting procedures could 
be applied. By “coherent” it is meant that all lines consist 
of three kinematic components (Table [3|) : two of them have 
similar z and are relatively narrow (FWHM~[95,210] and 
[470,560] km respectively). A third broad blueshifted 
component is moreover isolated in all lines, which is emit¬ 
ted by the outflowing gas. It has FWHM~[1400,1840 km 
s“^ and Vs ~[-490,-190] km s“^ . All the three kinematic 
components have line r atios consistent with type 2 AGN, as 
already pointed out bv IVillar Martin et al.l (|20I4l ). 

Other lines show clear evidence of a broad underlying 
component, although the triple Gaussian fit is not possible 
due to the low signal-to-noise and/or the complex blend with 


neighbour lines. We have calculated the ratio e = ^hroad 

-^narrow 

between the flux of the broad (outflowing) component and 
the rest of the line flux (which for simplicity we will name 
i^narrow, although it coutaius the narrow and intermediate 
components) for as many lines as possible. This gives a mea¬ 
surement of the relative contribution between the most tur¬ 
bulent outflowing gas and the more quiescent ambient gas. 

The ratio e is plotted against log(ncrzt), FWHM, IPiow 
and IPhigh in Fig. [8l A significant correlation is found with 
the critical density (rs=0.73 and p=0.003; panel A). No 
significant trend is found with the ionization potentials 
(rs=0.40 and p=0.I5 for e versus IPw and rs=0.35 and 
p=0.23 for e versus IPhigh)- 

An interesting case is the high ionization line 
[FeVH]A6087, for which a very broad component is isolated 
with FWHM=2460±340 km s“^ (Fig. [9} and Vs=-lb0±60 
km s“^. Thus, the line which shows simultaneously very high 
critical density and ionization potential, also shows the most 
extreme kinematics and the largest e = 0.96zb0.I0 of all opti¬ 
cal lines. At the other end, the lowest Ucrit (also low ioniza¬ 
tion) lines [NI]A5200 and [SII]AA67I6,6731, are the narrow¬ 
est (FWHM~300 km s“^) and have the smallest e ~0.1-0.2 
(Fig. El . 

The correlation between e and the FWHM suggests that 
the increasing broadening of the lines is due to the increas¬ 
ing outflow influence. The correlation with Ucrit shows that 
the outflow emission is relatively stronger in higher density 
gas. Ultimately, these results suggest that the FWHM versus 
Ticrit correlation is produced by the outflow in MRK 477. 

[Fell]Al.644 is the only forbidden NIR line for which 
the spectral decomposition could be applied. It consists of a 
narrow, spectrally unresolved component with FWHM<340 
km s“^ and a very broad and prominent underling com¬ 
ponent with FWHM=4770zb830 km (Fig. [TOl shifted 
by 3I5di23I km s“^ relative to the narrow core. It is the 
only line for which the broad component is not blueshifted. 
For comparison, Paa consists of two components with 
FWHM<260 km and 2240zb230 km s“^ respectively 
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Figure 9. Fit of the [FeVII] A6087 spectral profile. The presence of 
a very broad underlying wing (FWHM~2500 km s“^) shows that 
the high ionization Fe+® region also participates in the outflow. 
Flux in units of 10“erg s“^ cm“^ A“^. 
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Figure 10. Fit of the NIR [Fell] A1.644yL4m (top) and Paa (bot¬ 
tom) spectral profiles. Flux in units of 10“^® (top) and 10“^^ 
(bottom) erg s“^ cm“^ A“^. 


with the broad one blueshifted by -570zbll9 km s“^. As 
in Fig. O [FeII]A1.64/im is a clear outlier in all e diagrams, 
with e=1.3zb0.2, the highest of all lines. The line, therefore, 
has a dominant contribution of emission from the outflow¬ 
ing gas. The different behavior of this line will be discussed 
later (see below). 


3.4-2 Reddening and eleetron density 

(i) Reddening eorreetion. Reddening correction (c and E(B- 
V)) has been estimated using the Balmer ratios and 

for which we assume case B values 2.8, 0.47 and 0.26 
respectively, appropriate for gas densities n ~10^-10® cm“^ 
and e lectron temperatures Te ~10,000-20,000 K (Osterbrock 
Il989h . 

We show in Table [4] the E(B-V) and c values derived 
for the three kinematic components and the total line fluxes 
using the expressions: 


E{B - V) 
and 


1.99 X log 


(Ha/H/3)„, 1 

2.86 


[eq.l] 


r F 1 


r F 1 

Lh/jJ 

int 

iHhl 


[/(A)-/(H/3)] 


[eq.2] 


where “mt” and “ohs” denote the intrinsic and observed 
line ratios and /(A) — /(H/3) is give n by t he standard inter¬ 
stellar extinction curve (Osterbrock 1 19891 ) . The final values 
and errors take into account that negative c values are not al- 
lowe(0. Once c is determined, other line ratios are corrected 
for reddening using the equation (2). 

In spite of the complexity of the line profiles, the red¬ 
dening values derived from the three ratios are in reason¬ 
able good agreement. Erom Table 3 we conclude that the 
narrowest component shows little or no reddening while 
the intermediate component shows the highest reddening. 
Eor the broad component, the results are less clear, since 
the Balmer decrement ^ suggests higher reddening (E(B- 
V)=0.32zb0.06) than the and Hy ratios, which imply lit¬ 
tle or no reddening. We cannot discard problems with the 
Ha fit, which is severely blended with the [Nil] doublet. 
However, it might be a real effect. As we will see next, the 
broad component is expected to have high densities n ^8000 
cm“^, maybe up to ^10® cm“^, given the detection of very 
high critical density strong lines. Under these circumstances, 
the Balmer decrement can be enhanced due to collisional ex¬ 
citation of Hq . n > 5x10^ cm“^ will produce > 4 (e.g. 
Binette et al. Il993h . Thus, the inconsistency between the 
reddening values inferred for the broad component might 
actually indicate the existence of high densities in the out¬ 
flowing gas. 

We measure Pa(a/Br7=12.6zbl.4 from the NIR spec¬ 
trum. Taking errors into account, this implies no reddening 
or, at most, or E(B-V)^0.32. This is consistent with the 
range of values allowed by the optical decrements. 


This is the reason why in Table 3 the same line ratio with differ¬ 
ent errors (e.g. ^=0.23±0.04 and 0.23T0.06) can be associated 
with slightly different c and E(B-V) values. 
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and 3 in Keenan et al. 1996)). A trend is found such that the broader the lines, the higher the density. Thus, the outflowing gas has the 
highest density. 


(1) 

Comp. 

(2) 

Ha 

H/3 

(3) 

CHa 

(4) 

E(B-V)Ha 

(5) 

H7 

H/3 

(6) 

CH7 

(7) 

E(B-V)h7 

(8) 

as 

H/3 

(9) 

CHS 

(10) 

E(B-V)h5 

Narrow 

3.1T0.3 

0.14T0.12 

0.08T0.08 

0.47T0.03 

0.09T0.09 

0.04T0.03 

0.23T0.04 

0.37T0.37 

0.24T0.11 

Intermediate 

4.9T0.2 

0.70T0.04 

0.47T0.03 

0.38T0.06 

0.62T0.46 

0.42T0.19 

0.23T0.06 

0.40T0.40 

0.26T0.12 

Broad 

4.2T0.3 

0.49T0.08 

0.32T0.06 

0.57T0.09 

0.00 

0.00 

0.26T0.03 

0.15T0.15 

0.09T0.05 

Total 

4.0T0.2 

0.45T0.06 

0.29T0.05 

0.45T0.05 

0.21T0.21 

0.13T0.06 

0.24T0.02 

0.27T0.23 

0.17T0.08 


Table 4. Extinction correction (c and E(B-V)) determined from the Balmer ratios (columns 2, 5 and 8) for the three kinematic 
components and the integrated lines. 


(1) 

Comp. 

(2) 

(3) 

r2ohs 

(4) 

(5) 

(6) 

r3‘“* 

(7) 

n (cm ^) 

Narrow 

1.03T0.03 

0.05T0.01 

0.07T0.02 

0.007T0.005 

0.010T0.007 

-400-630 

Intermediate 

0.73T0.05 

0.17T0.02 

0.37T0.09 

0.08T0.01 

0.17T0.04 

-2000-4000 

Broad 

0.41T0.12 

0.24T0.05 

0.34T0.14 

0.08T0.07 

0.25T0.11 

>8000 

Total 

0.76T0.06 

0.13T0.01 

0.22T0.03 

0.05T0.01 

0.08T0.02 

-1600-3200 


Table 5. Density sensitive line ratios observed corrected (^’^^) for reddening, rl = [siljA673i ’ ~ [511^67:^^+6731 ’ ~ 

[Sll]^n6^673i • range of densities n has been calculated using rl. A trend is found such that the broader the component, the higher 
the density. 


(ii) Electron density n. We have constrained n with 
the ratio rl = [slrA673i ^Qsterbrock 1 19891 ) . The results are 
shown in Table [3 The apparently non perturbed ambient 
gas (the narrow component) has n ~400-630 cm“^. The in¬ 
termediate component, which shows intermediate kinematic 
and physical properties, has n ~2,000-4,000 cm Finally, 
the most extreme kinematic component, which is associated 
with the outflowing gas has the highest density n >8.000 
cm“^. 

We have further checked these results using r2 = 


[SII]4068 . o _ [sn]4076 fKeenan et al 

nlike rl, these ratios have a strong dependence on the elec¬ 
tron temperature Te and reddening. They are therefore less 
efhcient at constraining n, but they provide a useful test. 
The dependence of both ratios with n and Te can be seen in 
Fig. [11] 

To correct for reddening (eq. [2]), for each component 
we adopt c as the average of the maximum and minimum 
values obtained from the three Balmer ratios (Table [3|). The 
uncertainties are calculated using half the difference between 
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Figure 12. Several narrow Fe“^ emission lines are detected in the 
spectrum of MRK 477. [CaII]A7291 is also tentatively detected 
(third panel). 



Figure 13. Emission lines are detected from other refractory 
elements apart from Ca and Fe (Fig. 9), such as possibly Ni (top) 
and Mg (top and bottom). 
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Figure 14. Spectral fit, residuals and kinematic components of 
[OI] AA6300,6364. [SIIJA6312 (black small Gaussian) is marginally 
detected as a small excess on the red side of the [OIJA6300. The 
detection of [FeX]A6375 cannot be confirmed. [OIJA6364 is suc¬ 
cessfully fitted without the contribution of such a line. Color and 
line style code as in Fig. 5. 
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Figure 15. Top: “blue bump” near HeIIA4686 due to WR stars 
in the SDSS spectrum of MRK 477. Middle: the “red bump” 
(~5800 A) is also detected. Bottom: The optical spectrum of the 
WR galaxy SBS 1222+614 is shown for comparison (taken from 
Guseva, Izotov & Thuan 2000). Flux for this spectrum in units 
of 10“^® erg s“^ cm“^ A“^. 


these two values. The results are shown in Table [5] and Fig. 
ttn In this figure we show the n sensitive diagnostic diagrams 
log(rl ) ver sus log(r2) and log(rl) versus log(r3) (Keenan 
et al. 1 19961 ). The coloured squares mark the areas covered 
by the allowed range of extinction corrected line ratios for 
the different kinematic components and the integrated line 
fluxes, taking into account all the uncertainties. The figure is 
consistent with the results from rl. They support a trend for 
increasing density from the narrow to the broad component. 

In spite of the uncertainties inherent to fitting the mul¬ 


tiple kinematic components in all the lines involved in the 
determination of n and reddening, a coherent picture can 
be drawn: the ionized gas within a radius of 1.5 arcsec or 
~1.1 kpc (set by the size of the SDSS fibre) shows a gradient 
in physical and kinematic properties, which is apparent in 
the three kinematic components. The apparently non per¬ 
turbed ambient gas (the narrow component) has n ~400-630 
cm“^. The intermediate component, which shows interme¬ 
diate kinematic and physical properties, has n ~2000-4000 
cm“^. Finally, the broadest component, which is associated 
with the outflowing gas has the highest density>8000 cm“^. 
As mentioned above, the large ^ decrement might suggest 
densities as high as n >5x10^ cm“^ for the outflowing gas. 


3.5 Narrow optical Fe^ emission 

Broad Fe^ multiplet emission are promi nent features in the 
optical spectra of most type 1 AGNs. IPong et~M ] (l2010[) 
have demonstrated statistically that narrow optical Fe^ 
lines, either permitted or forbidden, are prevalent in type 
1 AGNs, but are completely absent in type 2 AGNs across 
a wide luminosity range, from S eyfert 2 galaxies to QS 02. 

MRK 477 is an exception. iHeckman et al.l (| 19971 ) de¬ 
tected [FeII]A8617 blended with some neighbour lines. We 
confirm the detection of more than 10 Fe^ emission lines 
(see Table[T]and Fig. [12]), besides several other features that 
might be blended or misidentified. The features for which 
the width could measured show FWHM~400-650 km s“^, 
suggesting similar kinematics as the NLR (Fig. (5] top pan¬ 
els). Other emission lines associated with highly refractory 
elements are also detected (Mg, Ni and, tentatively, Ga; Figs. 
[12] and [T3l) . 


3.6 Coronal lines 

Goronal lines are collisionally excited forbidden transitions 
emitted within low-lying levels of hi ghly i onized species 
(IPiow >100 eV; Rodriguez Ardila et al. 20061 ) . According to 
that exact definition, we cannot confirm the detection of any 
coronal feature in the optical spectrum of MRK 477, since 
the highest ionization lines we identify unambiguously are 
those produced by Fe^® (IPiow =99 eV). On the other hand, 
coronal emission is confirmed in the NIR spectrum with the 
detection of [SiVI]A1.963 (IPw =167 eV) and [SiX]A1.430 
(IPiow =351 eV) (Fig. 2] Table [2}. 

Althou gh [FeX]A6375 (IPi ow= 235 eV) was rep orted as 
detected bv IPe RobertisI (|l987l ) and IVeilleuxI (|l98 ^ , we can¬ 
not confirm this. Both groups used the same spectra in their 
analysis, one obtained in 1980 and another in 1985. The 1980 
spectrum shows a strong emission feature at the right A (see 
their Fig. 6, left), but this was discarded bv I Veilleu^ (| 19881 ) 
as an artifact. Although the 1985 spectrum shows a small 
excess near the [FeX] wavelength which the authors inter¬ 
pret as the detection of this line, our analysis shows that 
this is consistent with being due to the complex kinematic 
sub-structure of the [OI]A6374 line profile (Fig. [T^ . 

A broad faint and noisy feature is detected with cen¬ 
tral A =7888.8 A. Several lines are possibly contribut¬ 
ing, including MgIIA7877, [NiIII]A7889; [FeXI]A7892 and 
MgIIA7896 (Fig.[l3] top). Thus, the detection of [FeXI] emis¬ 
sion (IPiow=262 eV), which has been identified in some type 
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2 AGN (Penston et al. ll984 Rose et al. l2013l l cannot be con¬ 
firmed either. 

The highest optical ionization lines unambiguosly iden¬ 


tified are [NeV]AA3346,3426 (IPw =97 e V, ou t of the SDSS 
spectral range; reported by De Robertis Il987l i and the [Fe- 
VII] lines (IP w =99 eV) at A 5159, 6087 A (also reported 
by De Robertis 1987l b and lines at 3759, 4893, 5721 A. Also 
possibly [FeVII] at 5276 A, [ArX] at 5533 A, although other 
identifications cannot be discarded (Tabled]). 

With somewhat lower ionization level, [CaV] at 5309 
A (IPiow=67 eV) and several [FeVI] (IPiow=75 eV) lines at 
5146, 5176, 5335, 5485, 5631, 5677 A are also detected and 
possibly [FeV]A5726 (IPw=55 eV) and [FeVI] at 5426 A. 

[FeVII] A6086 is fainter relative to low ionization lines 
(e.g. [OI]A6300) than usually found in active galaxies 
with strong coronal emission, ^ 53^00 ^ =0.15zb0.01 for 
MRK 477. F or comparison, the Seyfert 1 and 2 galax¬ 
ies studied bv iRodrfguez Ardila et aP (|2006l i with detected 
[FeX] and [FeXI] lines show in general (4 out of 5 objects) 

‘ToilAetoo'' ™ the range ~0.6-5.7. 


3.7 The WR bumps 

I Heckman et al.l (|l997l i detected a broad emission complex 
around He HA4686, which is clearly appreciated in the SDSS 
spectrum. They fitted this so called “blue bump” satisfac¬ 
torily as a blend of the He H line together with other lines 
which are identified in Fig. d5](top). This unresolved bump 
is a blend of lines emitted mainly by late WN and early WC 
WR stars, although some contribut ion of early WN stars 
might be present (Schaerer V Vacca [l9^ i . iHeckman et "ahl 
lll997^ propose that it may be produced by an ensemble 
of about 30,000 WR stars (WN subtype), in which case 
MRK 477 would be a luminous (but not extraordinarily so) 
member of the class of WR galaxies. 

Another feature often identified in WR galaxies due 
to WR stars is the so called “red bump” due to broad 
CIVA5808 emission, emitted mainly by WCE stars. This is 
usually mu ch weaker than the “b lue bump” and was not 
detected bv iHeckman et al] (|l997l i . 

We show in Fig. [15] (middle panel) the MRK 477 spec¬ 
trum near the “red bump” and the optical spectrum of the 
WR galaxy SBS 1222+614 (Guseva, Izotov V Thuan 2000) 
for comparison (bottom panel). The appearance of the “red 
bump” in this and other WR galaxies is very similar to that 
in MRK 477. We thus believe that the “red bump” is de¬ 
tected in the SDSS spectrum of this QS02. Its very large 
width (FWHM ~80 A or ~4100 km s“^) rules out as its 
origin the broad wings due to the ionized outflow identified 
in many other emission lines. 


4 DISCUSSION 

We have analysed the optical and NIR spectra of MRK 477, 
th e nearest obscur e d qu asar. WR features first identified 
by I Heckman et al.l (|l997l i and confirmed here by the new 
detection of the “red bump” at ~5800 A (Sect. 3.7) show 
that it has undergone very recent star formation. 

The optical+NIR spectrum of MRK 477 is very rich, 
with ~100 detected emission lines (~90 in the optical). In 
spite of the lack of spatial information of the SDSS spectrum. 
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the spectral decomposition of numerous lines has allowed us 
to characterize, at least partially, the spatial structure and 
the gradients in the physical and kinematic properties of the 
gas. 


4.1 The origin of the FWHM versus Ucrit 
correlation 

As in many other type 1 and type 2 AGN, a significant cor¬ 
relation has been found between the FWHM of the lines 
and the critical density Ucrit (in log). Based on the analysis 
of e ight type 1 AGN spectra. Stern, Laor V Baskin (2014; 
IS14l hereafter) showed that this relation is consistent with 
the velocity field within the black hole gravitational sphere 
of influence (radius Rg), assuming the n oc R~^ relation 
implied by the equilibrium between the NLR clouds and 
the radiation pressure and that the emission of each forbid¬ 
den line is dominated by gas with n ~ ricrit- This explana¬ 
tion implies that beyond some threshold ricrit, n^rit, which 
depends on the luminosity in Eddington units rri, the line 
emission will be dominated by gas within and therefore 
is expected to show larger velocities than gas which kine¬ 
matics is dominated by the host galaxy. Eor MRK 477, the 
stellar veloc ity dispersion a* = 117 km s“^ (Zhang, Bian V 
Wang [ioosl i . Lboi and the values of Mbh quoted in Sect. 1 
imply Rg a nd 40 for log(MBH)=7.19 and 8.84 respec¬ 
tively (eq. 26 i n IS14li . This implies a threshold n°j.it = 10^'^ 
cm“^ (eq. 6 in lS14l) for log(MBH)=7.19. This value is much 
larger than the ricrit cm“^ where the EWHM 

starts to increase in MRK 477 (top-left panel in Eig. [5]), 
which makes the gravitational interpretation unlikely. Eor 
log(MBH)=8.84, = 10 ^'^ cm ^ is indeed consistent with 

the data. The success of this prediction supports the gravita¬ 
tional interpretation for the EWHM versus ricrit correlation, 
provided MRK 477 harbors such a massive black hole. On 
the other hand, the kinematic substructure of the lines raises 
doubts about whether the emission is truly emitted mostly 
by gas at n ~ ricrit- As an example, the [OHI] line flux is 
dominated by gas with n <4000 cm“^ (see Sect. 3.4). 

Alternatively, we propose that the outflow produces this 
correlation in MRK 477. Thisis supported by the correla¬ 
tions between e = ncvit (Sect. 3.3.1), which im- 

-^narrow ^ ^ 

ply that the outflow is relatively stronger and is responsible 
for the increasing line broadening at increasing densities. 
This correlation does not differ in any statistical sense from 
the usual EWHM versus ricrit one. The difference lies in the 
physical interpretation we propose. 

It is possible that the effects of ionized outflows explain 
the EWHM versus ricrit correlation in other AGN. 

Three kinematic components have been isolated in 
the main optical emission lines. The narrowest and in¬ 
termediate components, have similar z {Vs km s“^) 
and are relatively narrow (EWHM~ [95-210] and [470-560] 
km respectively). A third broad blueshifted component 
is moreover isolated in all lines with EWHM in the range 
~[1400-1840] km s“^ and Vs ~[-490-(-190)] km s“^. The 
narrow component is likely to trace the NLR ambient, non 
perturbed gas, while the broad component is emitted by the 
most turbulent, outflowing gas. The intermediate compo¬ 
nent traces gas of intermediate properties. All three compo¬ 
nents have line ratios consisten t with type 2 active galax - 
ies, as already pointed out bv IVillar Martin et al.l (|2014l i. 
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implying that they are spatially located within the quasar 
ionization cones. 

The difference in kinematics is associated with a dif¬ 
ference in physical properties. The broader the component, 
the higher the density. The sequence we hnd is n ~ (400-630) 
cm“^ for the ambient gas, n ~(2000-4000) cm“^ for the in¬ 
termediate component and n ^8000 cm“^ for the outflowing 
gas (see §3.3.2). Higher n of the outflowing gas is also sug¬ 
gested by the correlation between e = ^hroad ncvn for 

-^narrow 

different forbidden lines (see above). Density enhancement 
of the outflowing gas has b een fo und in some radio g alaxie s 
and QS02 (e.g. Holt et al. 1201 ll : Villar Martin et al. l2014ll . 


4.2 The spatial location of nuclear ionized outflow 

The isolation of the broad component in numerous emission 
lines demonstrates that the outflow involves gas covering a 
large range of ionization potentials and critical densities (at 
least IPhigh ~8 to 125 eV, log(ncrit) ~3.3 to 7.5). We try 
to constrain next the spatial location where it was triggered 
and how far its effects extend. 

Bennert et al. (2006a,b) found that the observed gas 
density of the NLR gas decreases with increasing dista nce 
to the AGN in Seyfert 1 and 2 galaxies (see also IS14h . A 
similar behavior is expected in MRK 477. The difference in 
kinematic and physical properties between the broad, inter¬ 
mediate and narrow components can be naturally explained 
if they are located at increasing distance from the AGN, with 
the outflow emission dominating at smaller distances, in the 
inner part o f the NLR or even closer (see also Villar-Martm 
et al. I 2 OI 4 I ). The outflow becomes weaker as it propagates 
outwards and reaches less dense regions, where lines of lower 
critical densities are preferentially emitted. Because the out¬ 
flow has lost some power, it drags less and less mass and its 
emission becomes weaker relative to the ambient gas, which 
dominates the line fluxes in the most distant, lower density 
regions. This can also explain the correlation of e with ricrit- 

If the rad iatio n pressure is responsible for the NLR den¬ 
sity gradient (|S14l) . the estimated n for the three kinematic 
components will help to constrain the distance from the cen¬ 
tral engine R at which each one is preferentially emitted. The 
expected behaviour of n with the ionizing luminosity Lion 
and R is predicted to be (|S14l b 

Tl — 7 X. 10 Lion,45 L5Q 

where Lion ,45 is the ionizing luminosity in units of 10^^ 
erg s“^ and R 50 is the distance in units of 50 pc. 

Typical ratios for luminous Syl and QSOl with 

-^ion 

Lbol in the range erg have median value of 

3.89 and standard deviation 4.36. Thus, for a ratio of 3.89, 
then Lion ,45 ~2.2 for MRK 477. Gonsidering the range of 
n inferred for each gaseous component, the broad (outflow¬ 
ing), intermediate and narrow kinematic components would 
be located at <220 (the upper limit is a consequence of the 
lower limit on n ^8000 cm“^), 375±65 and 880±100 pc re¬ 
spectively. This sets an upper limit on the distance at which 
the outflow has been originated of R <220 pc. If densities 
as high as ~5xl0^ cm“^ exist in the outflowing gas (Sect. 
3.4.2), then R <30 pc. For comparison, the radio jet extends 
at ~0.6 arcsec or 435 pc from the central engine. With all 


uncertainties involved, it is remarkable that these calcula¬ 
tions place the intermediate and narrow components near 
and beyond the edge of the radio jet respectively. This adds 
further support to the idea that the radio jet has originated 
the outflow. In this scenario, our results suggest that the 
outflow in MRK 477 is concentrated in the nuclear region 
and does not reach distances beyond ~fewxI00 pc. 

Alternatively, the higher density of the outflowing gas 
might be explained by the compression exerted on the out- 
flowin g gas by the radi o-jet i nduced shocks (Villar Martin 
et al. 119991 : Holt et al. l20Ilh . We note that this scenario 
is qualitatively different from the scenario where the NLR 
clouds are compressed by radiation pressure, as assumed in 
the models used above. The higher density would instead be 
located at or near the radio jet and decrease moving away 
from it. In this case a correlation would be expected between 
the morphology of the radio jet and the ionized gas (whose 
emission is enhanced due to shock excitation and/or density 
enhancem ent). T his is actually the case in MRK 477 (Heck¬ 
man et al. I 1997 I L Both the radio and [O HI] images show a 
bright knot or ridge of emission about 0.4 arcsec (~290 pc) 
to the northeast of the central source. This clearly shows 
that the radio jet is interacting with the NLR and enhanc¬ 
ing the emission at this location. Thus, an alternate scenario 
is that this detached knot is responsible for the bulk of emis¬ 
sion of the outflowing gas. If so, it must emit a wide range 
of lines, including [FeVH]A6087. High spatial resolution spa¬ 
tially extended optical spectroscopy would help discriminate 
between both scenarios. 

4.3 The origin of the NIR narrow [Fell] Al. 644 
emission 

We have found that this line is much broader than expected 
for its critical density and ionization potentials (Sect. 3.2). 
The underlying broad component has FWHM~4770±830 
km s“^, larger than any other line and contributes more 
than half of the total line flux (~57%). It also shows the 
largest e =I.3±0.2 and is the only line for which the broad 
component is not blueshifted. All these could be naturally 
explained in terms of reddening being much weaker in the 
NIR, so that both the approaching and receding parts of the 
expanding outflow are observed (unlike the optical lines). 
However, this is an unlikely explanation since Paa shows 
that the broad component is blueshifted by -570±II9 km 
s“^. 

Alternatively, it is possible that [Fell] Al.644 is relatively 
enhanced by the shocks induced by the outflow, more than 
other lines and at a less obscured spatial location (hence 
the non-blueshift). Outflow induced shocks have been often 
proposed as a relevant [Fell]Al.644 excitati on me chanism in 
active galaxies (e.g. Ramos Almeida et al. l2009l . Gontini et 
al. ?). As an example, ? studied the kinematics and exci¬ 
tation mechanism of the NIR emission lines in a sample of 
Seyfert I galaxies, including [Fell] A 1.257 and A 1.645 /xm. 
They found 3 out of 22 sources in which the [Fell] lines were 
the broadest, even broader than the coronal [SiX]AI.252 /xm. 
They suggest that in these objects [Fell] must arise from and 
additional source, partially formed in a region distinct from 
other low-ionization species, which they sugges t to be asso¬ 
ciated with shock excitation from the radio iet. lBlietz et al.l 
lll995h also suggest that the [Fell]Al.644 emission in the Sy2 
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galaxy NGC1068 is preferentially emitted by gas which is 
interacting with the nuclear outflow/jet and is possibly lo¬ 
cated at the interface between the outflow and the dense 
circumnuclear molecular clouds. 


4.4 Coronal lines 

Detection of coronal emission is confirmed in the NIR spec¬ 
trum, but not in the optical (Sect. 3.6). In the optical, the 
highest ionization lines are those emitted by Fe^®, which 
on the other hand are fainter relative to the low ionization 
lines than typically found in AGN with strong coronal emis¬ 
sion. Goronal lines are generally detected between just a 
fe w parsecs a nd a few hundred parsecs (e.g. Mazzalay et 
al. l2ninll2014h . They have been proposed to be formed in a 
region located at an intermediate distance between the clas¬ 
sical NLR a nd the broad-line region (BLR) (Muller Sanchez 
et al. l2006h . Alternatively, the coronal region might reside 
in th e inner wall of the dusty torus (Murayama & Taniguchi 
Il998l . Rose et al. l2013l j or in the low density surfa ce layer of 
the radiation pressure confined NLR clouds 1^. 

The detection of coronal emission in the NIR and its 
absence in the optical spectrum of MRK 477 suggests that 
this region is heavily reddened, partially hidden from our 
line of sight. Alternatively, since the coronal lines are usually 
highly nucleated (<1 kpc from the nucleus) and much more 
concentrate d than lower ionization NLR features (Rodrfguez 
Ardila et al. l201lh . it is also possible that they are diluted by 
the strong continuum contribution within the large aperture 
of the SDSS fibre (radius r\j 1 kpc). 

The [SiVI] and [SiX] lines hav e luminosities 40.0 and 
39.8 i n log and erg s~^respectively. [Rodriguez Ardila et al.l 
(|201l| j found that these lines display a narrow range in lu¬ 
minosity in Seyfert 1 and 2, with values for most objects 
located in the interval log(L)=39-40. MRK 477 is at the 
high end of this range. We measure [SiVI]/[SiX] = 1.6±0.1, 
also within the range measured for Seyfert 1 and 2 galaxies 
by those authors. 

4.5 Detection of narrow optical Fe^ emission 

We have identified more than 10 narrow optical Fe^ emission 
lines in the SDSS spectrum of MRK 477. The NIR [Fell] lines 
at 1.26 and 1.64 /am are routinely detected in type 2 S eyfert s 
(e.g. Ramos Almeida, Perez Garcfa V Acosta-Pulido [jpool j . 
However, to our knowledge MRK 477 is the first type 2 AGN 
with optical Fe^ line detections. 

To explain the absence of narrow optical F e^ in type 2 
AGN. while being prevalent in type 1 objects, IPong et al.l 
(|201ol j proposed that such emission is confined to a disk-like 
geometry in the innermost region of the NLR on physical 
scales of parsecs. This would be smaller than the obscuring 
torus and within the dust sublimation radius. Iron, which 
is a refractory element which easily condenses on to dust 
grains, is i n the gaseous phase in the absence of dust (Laor 
& Draine ll993l L In this scenario the narrow optical Fe^ 
emitting region is visible along our line of sight in type 1 
objects, but obscured by the extent of the dusty torus in 
type 2 counterparts. 

It is however difficult to picture a geometry such as 
this to explain the optical Fe^ emission in MRK 477. Broad 


permitted hydrogen and helium lines would be expected, 
unless a BLR did not exist, but this is not the case (Sect. 
1 ). 

We explore next whether the optical Fe^ features can 
be naturally explained by the intrinsic NLR emission. For 
this, we compare the measured integrated flux relative to 
H/3 of all Fe^ lines in the blue (4000-6000 A) and red (6000- 
7800 A) bands with the ratios predicted by photoioniza¬ 
tion models appropriate for the NLR conditions. We us e the 
radiation-pressure-confined NLR models described in IS14l . 
The assumption of these models implies that the line emis¬ 
sion is essentially independent of the ionization parameter 
at the slab surface Uq, as l ong as Uq 0.03 (see Fig. 3 
in S14; see also Dopita et al. l2002l and Groves et al. l2004| j . 
A continuous distribution of dusty gas as a function of the 
distance from the nucleus R is considered. This distribution 
is characterized by 77 , the power-law index of the dusty gas 
covering factor as a function of logarithmic unit of R. These 
types of models are successful in explaining various obser¬ 
vations of the NLR (see Dopita et al. 2002 and S14). 

We assume solar metallicity, an ionization sl ope o f 
<aion = —1.6 typical of luminous quasars (Telfer et al. I 2 OO 2 I L 
and ?7 = 0, i.e. a constant covering factor per log R, as sug¬ 
gested by the strong line ratios and the IR spectral energy 
distribution (S14). The photoioniz ation models are calcu¬ 
lated using CLOUDY (Ferland et al. l2013l L assuming hydro¬ 
static equilibrium. Depletion of refractory elements on to 
dust grains is taken into account, and the Fe^ ion is cal¬ 
culated according to the model described in Verner et al. 
(1999). 

The results for the FeII(4000 - 6000)/H/3) and 
FeII(6000 — 7800)/H/3) are shown in Table [6] considering 
both dusty and non-dusty models. The uncertainties on the 
observed ratios account for the errors on the flux measure¬ 
ments and the possible misidentihcation of some features 
(Table [T|). The models can reproduce successfully both the 
blue and the red bands. Therefore, the narrow optical Fe^ 
lines in MRK 477 can be naturally explained by the intrinsic 
NLR emission due to AGN photoionization. The high lumi¬ 
nosity and/or closeness of MRK 477 could be the reason why 
these lines have been detected, unlike other type 2 AGN. 

The effect of dust depletion produces ratios approxi¬ 
mately twice lower than those derived from dust-free models. 
Taking into account the uncertainties, the Fe^ band ratios 
do not discriminate between the two scenarios. Model pre¬ 
dictions of lines emitted by other refractory elements such as 
Mg, Ni, Ga (also detected for MRK 477, see Sect. 3.5) would 
help to infer whether the NLR contains dust-free gas. Only 
the [GaII]A7291 is predicted by cloudy. Galcium is much 
more sensitive to dust depletion than iron. The dust-free 
models predict [Ga II]A7291/H/3=0.08, while the line would 
be 1000 times fainter (and thus completely u ndetec table) 
in the dusty case (see also Villar Martin et al. ll966l L The 
[Ga II] line is tentatively detected with a ratio of 0.02±0.01 
relative to H/3 (Fig. [121 third panel). If this detection is con¬ 
firmed, it would imply that at least a fraction of the NLR 
gas is dust free. 
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(1) 

C FeIl4000-6000 ^ 

V H/3 ) 

(2) 

FeIl4000-6000 
\ H/3 ) 

(3) 

Modeli 

Dust 

(4) 

Model 2 

Dust-free 

(5) 

^ Felleooo-rsoo 

V H/3 ) 

(6) 

t FeIl6000-7800 
\ H/3 ) 

Modeli 

Dust 

ModeR 

Dust-free 

0.I4T0.03 

0.I5T0.I 

0.15 

0.36 

0.06T0.03 

0.05T0.04 

0.031 

0.053 


Table 6. Comparison between the measured integrated flux relative to H/3 of all Fe+ lines in the blue (4000-6000 A) and red (6000-7800 
A) bands with the ratios predicted by photoionization models appropriate for the NLR conditions (see the text). The superscripts 
and refer to the observed and reddening corrected ratios. Dusty and dust-free model predictions are shown. The narrow optical Fe+ 
lines in MRK 477 can be naturally explained by the intrinsic NLR emission due to AGN photoionization. 


5 CONCLUSIONS 


We perform a detailed spectral analysis and characterization 
of the NLR of the nearest obscured quasar (QS02) MRK 477 
at z=0.037 based on the SDSS optical and NIR H+K spectra 
obtained with WHT-LIRIS. 


(i) We confirm the h ybrid nature (AGN+starburst) pro¬ 
posed by Heckman et al. 19971 . based on the new detection of 
the “red bump” at ~5800 A due to WR stars, which implies 
that the system has undergone recent star formation. 

(ii) The optical spectrum of MRK 477 is rich in emis¬ 
sion lines, with ~90 detected features. In spite of the lack of 
spatial information of the SDSS spectrum, the spectral de¬ 
composition of numerous lines has allowed us to character¬ 
ize, at least partially, the spatial structure and the gradients 
in the physical and kinematic properties of the gas. 

(iii) Gas densities within the range ~1945l4gQ and up 
to >10^ cm“^ are confirmed in the NLR of MRK 477. 


(iv) As in many other active galaxies (AGN), a signifi¬ 
cant correlation is found between the lines FWHM and the 
critical density log(ncrit)- We propose that this correlation is 
caused by the outflow and its impact on the gas kinematics. 
This could be the case in other AGNs. 


(v) MRK 477 is an example of a radio quiet powerful 
AGN where negative feedback (the nuclear outflow) can be 
dominated by the radio structures. The outflow emission has 
been isolated in many emission lines covering a large range of 
ionization potentials and critical densities (from [OI]A6300 
to [FeVII]A6086). The outflowing gas, which is concentrated 
within R^fewxlOO pc from the central engine, is >13 times 
denser {n >8,000 cm“^) than the ambient non perturbed 
gas (n ~400-630 cm“^). Marginal evidence is found for den¬ 
sities as high as n >5x10^ cm“^ in the outflowing gas. It is 
possible that the density enhancement is due to the gas com¬ 
pression produced by jet induced shocks. This is supported 
by the correlation between the rad i o and [OIII]A5007 mor¬ 
phologies found bv iHeckman et all (|l997l ). 

Alternatively, the density enhancement is not related to 
the jet. Instead, it might be a reflection of the NLR intrinsic 
density gradient, consequence of the gas being compressed 
by radiation pressure. In this scenario, and based on the 
comparative study between the density and the kinematic 
properties of the outflowing and the ambient gas, we con¬ 
clude that the outflow has been generated at <220 pc (pos¬ 
sibly at <30 pc) from the AGN. We find evidence of how 
its effects weaken as it propagates outwards, following the 
NLR density gradient. Beyond the radio jet edge, the gas 
emission is dominated by ambient less dense, non perturbed 
gas. This adds further support to the idea that the radio jet 
has triggered the outflow. 


(vi) The [FeII]A1.644/im line presents a very different 
behaviour than the rest of the emission lines. It shows the 
most extreme effects of the outflow, with an underlying 
broad component of FWHM=4770zb830 km s“^. Its proper¬ 
ties suggest that its emission is enhanced by shocks induced 
by the nuclear outflow/jet and is preferentially emitted at a 
different, less reddened spatial location, maybe the interface 
between the outflow and the dense circumnuclear molecular 
clouds. 

(vii) More than 10 narrow optical Fe< emission lines 
have been detected in the SDSS spectrum of MRK 477. To 
our knowledge, this is the first type 2 AGN with such a de¬ 
tection. We show that these lines can be explained as the 
natural emission from NLR gas photoionized by the AGN. 
Emission lines associated with other highly refractory ele¬ 
ments (Mg, possibly Ni) are also detected. If the tentative 
detection of the [Ga II]A7291 line is confirmed, this would 
imply that at least part of the NLR gas is dust-free. 

(viii) Goronal line emission is confirmed in the NIR, 
but not in the optical SDSS spectrum. The coronal region 
might be heavily reddened, partially hidden from our line 
of sight. Alternatively its optical emission might be diluted 
due to the large SDSS fibre aperture. The coronal region 
also participates in the outflow. 

(ix) Paa is spatially extended along the radio and [GUI] 
emission axes, up to a maximum radial extension of ~ 1.5 kpc 
from the AGN. 
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